The capture into parametric autoresonance* 

O.M.Kiselev, t S.G. Glebov* 



Abstract 

In this work we show that the capture into parametric resonance 
may be explained as the pitchfork bifurcation in the primary para- 
metric resonance equation. We prove that the solution close to the 
moment of the capture is described by the Painleve-2 equation. We 
obtain the connection formulas for the asymptotic solution of the pri- 
mary parametric resonance equation before and after the capture using 
the matching of the asymptotic expansions. 

Introduction 

We study the primary parametric resonance equation: 

*e0'+(a(fl) + |0| 2 )0-0* = O. (1) 

Our goal is to develop an asymptotic theory of the capture into the para- 
metric autoresonance with 0(6) = —9. We give the asymptotic analysis of 
the capture into the parametric autoresonance and connect the asymptotic 
formulas for the solution of before and after the capture. 

The capture in a nonlinear (not parametric) resonance was studied in 
0. The capture in the nonlinear resonance associated with loss of the sta- 
bility and slow crossing through the separatrix (see PUIS])- Autoresonance 
phenomenon was supposed for accelerators of relativistic particles in [H E] . 
Later the autoresonance was founded in many different oscillatory and wave 
processes (see review jB]). Mathematical approach to the capture into the 
autoresonance was considered in [7] . The capture into autoresonance accom- 
panies with an bifurcation or separatrix crossing in the primary resonance 
equation (see |Hj). 
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It is well-known the parametric resonance leads to the exponential grows 
for the solutions of the linear equation (see [§]). However the nonlinear 
terms lead to the unbalancing of the frequency of the external force and the 
frequency of the oscillations for the solution of the nonlinear equation. It 
changes the amplitude of the oscillations for the value of the order of square 
root of the perturbations. This phenomenon was shown in ^U] by analysis 
of the solution for primary parametric resonance equation in the form: 

R' + R sin(2z/) =0, v - a - R 2 + cos(2z/) =0, (2) 

where a = const. 

Later equation (J2J) was studied in ^] EI], when a ^ const. It was 
shown that the autoresonance phenomenon takes place for the primary para- 
metric resonance equation when a < and the modulo of the coefficient \a\ 
grows with respect to 9. 

We consider the primary parametric resonance equation in the form (Q). 
One can obtain this equation after the substitution: = Rexp(ii>). There 
are two reasons to investigate equation with a(8) = —6. On the one hand 
the case cr(6) = —9 is the simplest one and it contains the autoresonance 
phenomenon for parametrically driven systems. On the other hand this case 
saves most of essential features of the solution for equation (£[]). 

Our analysis gives opportunity for an another point of view on the sub- 
ject. The capture into the parametric autoresonance may be considered as a 
loss of stability in the pitchfork bifurcation. The dynamic theory of the pitch- 
fork bifurcation for ordinary second order differential equations with slowly 
varying coefficient was considered in [TH UK] ITK] . The bifurcation when two 
centers and saddle coalesce is called by supercritical pitchfork. The solution 
in bifurcation layer approximates by solutions of the Painleve-2 equation (see 
|14j). Later the connection formulas was obtained for the solution before and 
after the supercritical pitchfork bifurcation for a asymptotic solution of the 
special form of perturbed Painleve-2 equation with dissipation and slowly 
varying bifurcation parameter (see [TH]). 

The supercritical pitchfork bifurcation for equation 

A' = i(XA - \A\ 2 - 5A*) 

was considered in . Here A and 5 are two parameters. We study the same 
equation but the parameter A varies slowly. 

In this work we show that the capture into parametric resonance may be 
explained as the pitchfork bifurcation in the primary parametric resonance 
equation. We prove that the solution close to the moment of the capture is 
described by the Painleve-2 equation and obtain the connection formulas for 
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the solution of the primary parametric resonance equation before and after 
the capture using the matching of the asymptotic expansions. 

The contents of the paper is as follows. The first section contains the 
discussions based on the numerical and qualitative analysis of the problem. 
This section contains the short description of the main result also. The WKB- 
asymptotic expansion for the solution of (JH) before the caption is obtained 
in section 2. In section 3 we match the WKB- asymptotic expansion and 
the asymptotic expansion based on the increasing solution of the Painleve-2 
equation. This equation is used into a small resonant layer. Our analysis 
of the solution for perturbed Painleve-2 equation is based on the connection 
formulas for the Painleve-2 transcendent (see (T^IES])- The WKB-asymptotic 
expansion for the captured solution of (|TJ) is constructed and matched with 
the increasing solution of the Painleve-2 equation in section 4. 

1 Discussion of the problem and main result 

In this section we discuss the capture into the parametric resonance for the 
solution of (0) and present the numeric and qualitative analysis of the cap- 
ture. We formulate the main result of the work at the end of this section. 

1.1 Numeric analysis 

Let us consider two numeric solutions for equation (Q) (see figure These 
solutions differ at the initial moment only. The first solution (left) corre- 
sponds to the initial moment at t = — 2 and the second one (right) corre- 
sponds to t = —2.01. We see the solution of the equation is very sensitive 
with respect to an initial data. 

The picture of |0| 2 for the oscillations looks like the figure |U for both 
solutions. 

Figure El shows the moment of the capture into the autoresonance. This 
moment is closed to t — — 1. 

The illustrated phenomenon is the scattering on the parametric autoreso- 
nance. Our goal is to obtain the solution for the scattering problem in terms 
of the asymptotic theory. 

1.2 Qualitative analysis 

Here we present the qualitative analysis for the equation with the "frozen" 
coefficient a = — T: 




(3) 
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Figure 1: the figure shows two solution of (JTJ). Left curve shows the solution of 
the Cauchy problem for the equation (JTJ) at t = —2.01, 3ft(0) = 0.02, 3(0) = 
0, e = 0.01. Right curve shows the solution of the Cauchy problem for the 
equation (0) at t = -2, K(0) = 0.02, 3(0) = 0, e = 0.01. 

The trajectory of the solution for the equation with the varying coefficient 
a(9) = —9 is locally closed to the trajectories of the equation (JHl) at T = 9. 
Therefore the solution of the equation with the frozen coefficient gives a local 
behavior for the solution of 

The Hamiltonian for the equation is: 

<n = -\w + m 2 + \m 2 + 2 )- 

It easy to see that there exists only one center when T < — 1. There exist two 
centers and one saddle when — 1 < T < 1 and three centers and two saddles 
at T > 1. We show the phase portraits for the equation with different value 
of the parameter T on the following figure. 

These three pictures give the conjecture about the numeric solution that 
was shown on the figures Q and El On figure Q one can see the non resonant 
solution when t < —1, the capture when —1.1 < t < —0.8 and the captured 
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-2 -1.5 -1 -0.5 

Figure 2: This figure shows |0| 2 for the solutions of the Cauchy problem for 
equation JTJ) at t = -2, 3fc(0) = 0.02, 3(0) = 0, e = 0.01 and at t = -2.01, 
9fc(0) = 0.02, 3(0) = 0, £ = 0.01. 




Figure 3: On the left picture T j -1, on the middle picture T=0, and on the 
right picture T i 1. 

solution when — 1 < t. Before the capture the solution oscillates around the 
unique center for the T < — 1. When T > — 1 we see the captured solution 
that oscillates around one of two centers from the second picture in figure El 
The left trajectory is captured by the left center from the second picture on 
figure El and the right trajectory on figure His captured by the right center 
on the second picture of figure El The bifurcation at the T = 1 does not 
effect on the captured solutions. These solutions remain close to the same 
centers. Below we obtain asymptotic solutions that illustrate the qualitative 
and numeric analysis. More over we obtain the formula that connects the 
asymptotic solutions before the capture and the solutions captured by the 
left or the right center. The formula is obtained by the matching method 

129- 
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1.3 Main result 

Let us formulate the main result. 

Theorem 1 Let the asymptotic solution of the primary resonance equation 
be 



9-1 

e + T 



sin 



9-1 



cos 



W 



as 9 < — 1, 



where 



and 



UJ 2 
S = ~ + ^1,0 + "1,0 



29 + 2 In 



9-1 



9 + 1 



u = -e^ffT^l - - ln(0 + VW^l). 



The constants ck^o and </? 1)0 are the parameters of the solution and 

Vi,o 7^ ^«i,o ln ( 2 ) - j ~ ar S f r (* a i,o/ 2 )J + X7r(mod(27r)), x = 0, 1, 



£nen £ne asymptotic solution has a form: 
$V\S t 6,e) ~ (-l) j 
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~ v/TT^ + £ 1/2 A ,o (77^= cos(5) + i^I+^sin^) ] 
<A>,o + Al fl (h + |e a + i(112 - 89)VT+9 + ? ln(l + 6)j 



in the domain 9 > — 1. ifere 

j = 2, 3, the constants <£>o,o anc ^ A),o are the parameters of the asymptotic 
solution. The parameters are defined by 



0,0 



In 



7T 



3|3(p)| 



7T 7 



V^o.o 



4-arg T(< ) -arg(l+p : 



where 



p= y / exp(vra? i0 ) - 1 exp fi-a£ ln(2) - -iarg (r(ia? (0 /2)) -2^1,0 
In formula j = 2 as ^(p) > and j = 3 as < 0. 
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2 WKB-solution before the capture 



Let us construct the solution of the WKB-type before the capture. The qual- 
itative analysis shows that the capture into the parametric resonance occurs 
at 9 = — 1. When 9 < — 1 the solution of equation (J3J) has an unique center 
at ip — 0. This fact prompts that equation has the oscillating solutions 
and they are closed to ip = 0. We construct the asymptotic expansion for 
the solutions of such type in this section. 

2.1 The WKB solution closed to zero 

Let us consider the WKB-solution in the form: 



oc 

-1/2 



^V-VnM). (5) 



(p = e 

n=l 

The leading-order term of the asymptotic expansion is the solution of the 
equation: 

lu'd^x - 00i - <p{ = 0. (6) 
The higher-order terms of (jSJ) are solutions of the equation: 

iu(j)' n - 9(j) n - (f)* n = f n , n = 2, 3, . . . , (7) 

where 

f n = -id 6 (f)n-l - idetpn-xdgCpn + h n . 

Here h n is a polynomial of the third order with respect to <f) m , m < n. 



The basis of the solutions for the linear equation To solve equation 
(JBJ) and the equations for the higher-order terms we construct the basis of 
two linear independent solutions for equation (JHJ) as the WKB-approximation. 
Define the basis solutions for equation (jUJ) by V\ and i^. Let the first solution 
be 

vi(t, 9) = cu' cos (r) — i{9 + 1) sin (r) , 

where, r = u/e, 

(cu') 2 = 9 2 -l. 

The Wronskian of two linear independent solutions of equation (jUJ) is con- 
stant. Therefore we construct the second solution as: 

— 1 i 

V2 ^ 9> > = 7e + T) sin M ~ 17' cos ( r ) ' 
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The Wronskian for these solution is: 

w(vi, v 2 ) = ViV% — v{v 2 = 2i. 
The formula for the general solution of (JHJ) has a form: 

v(t,6) = C 1 (9)v 1 {6,S) + C 2 {6)v 2 (t,9). 
Using the WKB approximation we specify the C\{9) and C 2 (9) such that: 



C 2 = C 2>oy /(-0 - 1) 



Here C1.0 and 6*2,0 are arbitrary constants. 
We use the functions 



and 



Wl(T, 



w 2 (t,9) 



9-1 
9 - 1 



COS T 



sin it) 



w. 

9 + 1 
9-1 

9 + 1 



sin r 



cos r 



I J \9-l 
as the basis of the linear independent solutions of equation (JHJ). 



2.2 Constructing of the WKB solution for nonlinear 
equation 

The main term of the WKB asymptotic solution for (0) we construct in the 
form: 



4>i(s, 9) = a 



0-lV 



sins 



9 



COS [S 



9 + 1) w 

where s = uo/e + ip and arbitrary constants a and are parameters of the 
solution. These parameters are represented by asymptotic series: 

00 00 

k=l k=l 

It is easy to see that the solution of (JZJ) is bounded if the right-hand side 
is orthogonal for two linear independent solutions of the linearized equation. 

f n (z, 9)wt(z, 9) + /*(«, 9) Wl (z, 9)dz = 0, 



(9) 



2- 



f n (z, 9)w*(z, 9) + f*(z, 9)w 2 (z, 9)dz = 0. 



Let us construct the bounded solution of (JJJ) at n — 2. The right-hand 
side is 

fi = _ l^ilVi = -ideaxd a <j)i - iaid s (pi - |0i| 2 0i. 

Formulas © give the equations for the functions a±(9) and (fi(9): 

cWi = 0, 2<9 e (/?i + a 2 |0i| 2 = O. 
Using the formula for (pi we obtain: 

J6 + l) 2 + (6- l) 2 

2 ^' +a ' (s+V-D =0 - 

The solution of the equation has the form: 



Analogously, the next terms of the WKB-asymptotic expansion are con- 
structed. The equation for the n-th correction term of the WKB-asymptotic 
solution has a form: 

%J$ n - e<p n -<p* n = f n , (10) 

where 

/„ = d e a n -id a (j)i + ad e (pn~id s (t)i + K. 

Here h n depends on the lower order terms of the WKB asymptotic expansion. 
These equations allow to determine coefficients a n and ip n . 

Asymptotic behavior close to the turning point The WKB asymp- 
totic expansion is not valid at the turn point 6 — — 1. As 6* ^ — 1 — the 
asymptotic behavior of the coefficients for asymptotic series (0) are repre- 
sented by following formulas. For the leading-order term we obtain 

0i ~ « 1)0 (j^l) ' sin (^(l-#) 3/2 -Ko) 2 m(l+0)+^ 1)O ) , - -1-0 

where a.\$ = const and ipi^ = const are the parameters of the solution. 
The second term of asymptotic series (0) has a behavior: 

2 = o((-i - ey 1 ^) e->-i-o. 

The recurrent calculations give 

n+1 = O(0 n (-1 - #r 3/2 )), - -1 - 0. 
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The domain of validity Expansion (jSJ) saves the asymptotic property 
with respect to e when 

e 4>n+l 1 a in 

— - — <1, as — > — 1 — 0. 

It yields the domain of validity for ((HJ) 

e~ 2/3 (-l-6) > 1. 

3 The Painleve layer 

To construct the asymptotic solution when 9 is closed to —1 we use the scaled 
variables: 

8 + 1 = e 2/3 r], (f) = e 1/3 x(r], e) + is 2/3 y(r], e). 

This change of the variables leads to the system of the equations for real and 
imaginary parts of the function <fi: 

x' + 2y = e 2/3 {r] - x 2 )y - e 4/3 y 3 , 

y + (77 - x 2 )x = e 2/3 y 2 x. (11) 

3.1 The asymptotic expansion in the Painleve layer 

We construct the solution of system (|11|) of the form: 

00 00 
x(r ] ,e) = Y,e 2n/3 x n (r ] ), y(r), e) = £ e 2n / 3 y n ( V ). (12) 

n=0 n=0 

The leading-order terms of (j!2)l are solutions of the system of equations: 

x'o + 2 Vo = 0, 
y' +( V -x 2 )x = 0. (13) 

The higher-order terms are determined from the system: 

x' n + 2y n = H£\ 
y' n +( V -3x 2 )x n = H^. (14) 

Here functions Hn^ and H^f 1 depend on 77 and the lower-order terms of 
asymptotic expansion f|12j) . For example 

h[ 1] = (77 - x 2 )y , H[ 2) = ylx ; 
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#2 = (V ~ x l)Vi - ^x Q y Q Xi - yl, H%> = 2y y 1 x + y%Xi, 

The leading-order term of the asymptotic expansion xq is a solution of 
the second order equation: 

x'q + 2(— 77 + x 2 )x = 0. (15) 

The higher-order terms of f!12l) satisfy to the linear differential equation 
of the second order: 

x" n + 2C-77 + 3xg)z n = ^, h n = + 2^#«. (16) 
The right hand side of (JTfij) has the following structure: 

n— 1 n— 1 

^ n = # f + 2^ + 2 Y, d x MH + 2 £ d vM ] y' k . 

k=0 k=0 

The substitution x' k = —2y k + and y' k = (3xq — r])xk + ^ 2 ' ) gives: 

n-l 

h n = + 2 ^^) + 2Y J d x M\-2y k + /#>) + 

fc=0 

n-l 

2E^^ 1) (( 3l o-^ + ^ ) )- ( 17 ) 

fc=0 

The substitutions 

7] = 2 1/3 z, x (r}) = 2 1/3 iu(z) (18) 

reduce (jT5j) to the Painleve-2 equation: 

u" - - 2u 3 = 0. (19) 

Using the substitution x n (rj) = 2 1 ^ 3 iu n (z) we obtain the perturbed linearized 
Painleve-2 equation: 

v!' n - zu n - 6u 2 u n = H n . (20) 

The Painleve-2 equation is integrable by the isomonodromic deformation 
method [21]. The solution u(z; a, 0) of the equation is called by the Painleve 
transcendent that depends on two parameters a and (p. It is known that the 
real solution of (|T9*|) does not have the singularities as z G K Therefore 
the main term of the asymptotic expansion is presented by the Painleve 
transcendent u(z,oc,0). 
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The homogeneous linearized Painleve-2 equation: 

w" — zw + 6u 2 w = 
has two linear independent solutions: 

u\ — da.u(z, a, 0), u 2 = dq,u{z,oi,(p). 

These solutions allow to present the higher-order terms of the expansion (|12|) 
in the form: 

x n = B^XM + B^X 2 (z) + Xl {z) f MlMl dC _ 

J no vv 

XM rwm dc (21) 

Here = const 7^ is the Wronskian of the solutions Xi = iui(2~ 1 ^ 3 r]) and 
X 2 = iu 2 (2~ 1 ^ 3 rj) for the linearized equation, Bn\ B$ and r] are arbitrary 
constants. 

The asymptotic behavior on the left-hand side of the validity inter- 
val In this paragraph we determine the domain of validity for asymptotic 
expansion (|T2*|) and match the expansion with the WKB-asymptotic expan- 
sion that was obtained in section |21 Therefore we should determine the 
asymptotic behavior of the Painleve-2 transcendent u(z, a, (p) as z — > —00. 

The asymptotic expansion of the Painleve-2 solution as z — > —00 has a 
form (HUES]): 

u{z) = ta(-z)- 1 / 4 S m(^(-z) 3/2 + ^a 2 ln(-z) + <p) + o((-2)" 1/4 ). (22) 

Here a and (p are parameters of the solution of Painleve-2 equation. 

Formula (J22J) give the asymptotic behavior for the solutions of the lin- 
earized equation U\ and U 2 as z — > —00: 

U x = d & u ~ ^Q,2i n (_^(_ z )-V4 C0S (^_^3/2 + 3 a 2 \n(-z) +0), 

U 2 = d$u ~ ia(-z)-^ C os(^(-^) 3/2 + ^a 2 ln(-z) + 0). 

Using these asymptotic behaviors we obtain the asymptotic formulas for 
the higher-order terms as z — > —00: 

x n = A-XM + B~X 2 ( V ) + X-(z). (23) 

Here A~ and B~ are constants. They are obtained by matching procedure 
with external asymptotic solution below and X~ (z) is the particular solution 
of the equation for the n-th term of the asymptotic expansion x n . 
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The left border of the interval of the validity for the Painleve layer 

Using the asymptotic behavior of Xq as rj — > — oo and the asymptotic behavior 
of X\ and X2 we obtain the asymptotic behavior of the higher-order terms. 
The equation for x\ is 

x'[ + {-ri + x 2 ) Xl = 0({-7]) 2 - 1 / 4 ), 7] -> -00. 

The particular solution has an order: 

xr = o((-^) 3/4 ), v ^-oo. 

The higher-order term satisfies the equation of the form 

X n + (-V + x l) x n = 0((-?7) 2 X n _i). 

By the sequential calculations we obtain 

X; = 0((-77f- 1/4 ). 

This estimated value for the higher-order terms give the left border for 
the interval of validity for (|12|) : 

p 2(n+l)/3„ r 2 / 3 X"" 

or 

-e 2/3 r] < 1, 77 < 0, e -> 0. 



3.2 Matching with the WKB asymptotic expansion 

The domains of validity of asymptotic expansions © and (fT!?J) are intersected. 
The matching with the WKB asymptotic expansion from section 14.21 allows 
to obtain the parameters a and (p for the solution of the Painleve-2 equation, 
for example: 

a = 0:1,0, <P = <Pi,o- (24) 



Asymptotic behavior on the right-hand side The theorem by Its- 
Kapaev-Belogrudov [T%| ITH] gives the asymptotic behavior of solution (122(1 
for the Painleve-2 equation. 

Theorem 2 (by Its-Kapaev-Belogrudov) If a solution of the Painleve 2 
equation has asymptotic behavior \2ty) as z — > —00 and 



£1 = -a ln(2) - - - arg ( T( — ) ) + K7r(mod(27r), k = 0,1, 
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oo, 



then the solution of the Painleve-2 equation is: 

<z) = ^Lz- 1 '* exp (-^ 3/2 ) (l + o(l)), z 

where a 2 = exp(7ra 2 ) — 1, and sgn(a) = 1 — k/2; 
Otherwise the absolute value of the solution increases: 

u{z) = iz^l ± z(2z)^pcos (^ z 3/2 - ^p 2 \n(z) + v^j + o(^" 1/4 ), (25) 

where p > and < v < 2tt . The parameters p and v is defined by a and 
<p: 

1, fl + \p\ 2 ' 



p 2 = -\n 



ir \3\5s(p)\ 



71 7 

v = - 

where 



v = —j + ln(2) - arg(r(ip 2 ) - arg(l + p A 



p = (exp(vra 2 - 1) 1/2 exp (i^oi 2 ln(2) - - i arg(r(z'y)) - itp\ , 

and the sign "+" is defined by 5s (p) < 0; 
if p = 0, then 

fz 2~ 5 / 2 

u{z) = ±i^-±ijj= + 0{z-W). (26) 

These formulas give the asymptotic solution of the system of equations (|11|) as 
rj — > oo. The solution of the scattering problem for the system of equations 
(jll|) allows to solve the problem of the capture into the resonance for the 
small amplitude solution of the primary resonance equation. 



The right border of the interval of the validity for the asymptotic 
expansions in the Painleve layer The formulas (|2*H|) . (|2T)|) and substi- 
tutions (fTHjl give the leading-order term of the captured into the resonance 
asymptotic solution. The solution has the behavior: 

Xo{rj) ~ =F2 1/3 v /r/, 77 -> oo. 

The nonlinear terms in the perturbed Painleve-2 equation (jllj) are 

x n ( V ) = o((- v r +i / 2 ). 

Then asymptotic expansion (|T2*j) is valid in the domain 

e 2/3 rj < 1. 
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4 The captured WKB-asymptotic solution 



The previous section shows the asymptotic solution of equation (JTJ) increases 
as ±Vd. When 9 — > —1 + the leading-order term of the asymptotic ex- 
pansion in the Painleve layer corresponds to the one of the centers for the 
equation with frozen coefficient. Equation (Q) has two slowly varying solu- 
tions as — 1 < 9 < 1. In this section we construct the formal asymptotic 
expansion for these solutions and for slowly varying solutions as 9 > 1 also. 



4.1 Slowly varying solutions 

Let us to construct the solution of the equation (PJ as follows: 

oo 

U(9,e)= y *Te n U n (9) (27) 

n=0 

The main term of the asymptotic expansion is defined by an algebraic 
equation: 

-9U + \U \ 2 U -U* = 0. (28) 
Using the complex conjugated function we obtain: 

[U o + U o *][U o -U *]=0 

This formula shows that the leading-order term of the asymptotic expansion 
is pure real or pure imaginary one. The real terms are 

f/ (1) = 0, f/f 3) = ±vT+~0, (29) 

and the imaginary terms are: 



[/ (4 ' 5) = ±iVW^l. (30) 

The algebraic equations for the higher-order terms Un\ j = 2,3,4,5 are 
follows: 

[-9 + 2\U^)U^ + [(U^) 2 - 1](U^Y = Q n . 

Here 

k+l+m=n 

where Cuim is a constant. 

Let us represent the linear equation as a system of the equations for the 
real and imaginary parts. The determinant of the system of equations is: 

A&> = 4(1 + 0), j = 2, 3 
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and 

A {j) =4(1-9), j = 4, 5. 

Then the system for the higher-order terms is solvable if j = 2, 3 and 9 > — 1. 
Otherwise if j = 4, 5, then the system for higher-order terms is solvable when 
9 > 1. 

When 6* — > oo the order of n-th term of asymptotic is 0(9~ l l 2 ), as n = 
1,2, ... . Therefore the asymptotic expansion (}2Tj) is uniform with respect to 
9 when 9 > — 1 for j = 2, 3 and when 9 > 1 for j = 4, 5. 



4.2 WKB-asymptotic expansion close to the slowly vary- 
ing centers 

Here we construct the WKB-asymptotic expansion that oscillates close to the 
slowly varying centers The WKB asymptotic expansion has a following 
form: 

oo 

m e) = V® (9, e) + Y, e k ' 2 ^\S, 9, e). (31) 

fc=i 

Here 

oo 

S = Q(9)/e + J2^ k/2 M0), 

k=0 

The equations for the higher-order terms are: 

itfds&P + [-9 + 2|C/ (i) | 2 ]^' ) + [(U^) 2 - 1]$?* = G n , (32) 
G n = -i ds^^-ide^- C^U®)*®*®*®*, 

l+k=n—2 k+l+m=n 

where CkimiU^) are the polynomials. 



WKB-solutions for the equation in the variations The homogeneous 
linearized equation of (}3*2*j) is called by equation in the variations: 



in'd 01 V + (-9 + 2\U {j) \ 2 )V + [U U)2 - 1}V* = 0. 

Two linear independent solutions of this equation have the following forms. 
The first asymptotic solution is: 

V X {9 X , 9) = Q' cos(0i) +i((-9 + 2\U ij) \ 2 ) + [U ij)2 - 1] ) sin(^). 
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Here 9 1 = fi/e, 

ft' = >/ [ — 6> + 2|C/0')| 2 ] 2 — I C/0'> 2 — 1 1 2 , 
The second asymptotic solution is: 

V ^ 6) = [-6 + 2\UV)\*] + WW-l) Mdl) " 

)|2 

cos(^i). 



([/(i)2 _ !) + (_6> + 2|t/Ci 



Q'([-9 + 2\U^\ 2 ] +$l(U^ 2 - 1)) 
The Wronskian of these solutions is: 

w(v u v 2 ) = v x v; - v*v 2 = 2i 

Using explicit formulas for the asymptotic expansion of U^forj = 2,3 
as e — > one can obtain: 

ft' ~ 2v / rT^, 

7i ~ 2Vl + # cos(0i) + 2%{9 + 1) sin(0i) 

V2 ~ — ; sinf^i) — % — ; cos(^i). 

" 2(1 + 0) V 1J 2v^TI 

The first correction term of the WKB asymptotic expansion Let 

the first correction term of the WKB-asymptotic solution <&i\s,9,e) be 

$?(S,9,e)=A(9,e)V 1 (S,9), 

where 

00 

A(9,e) = J2^ k/2 Md)- 

k=0 

The first correction term of the asymptotic expansion ()31|) has two series 
of parameters A k (9) and & G N. These parameters are defined by 

usual way in the WKB theory. The parameters are solutions of the anti 
resonant equations: 

2W G k (z, 9)V 1 *(z, 9) + G* k (z, 9)V 1 {z, 9)dz = 0, 

(33) 

G k {z, 9)V*{z, 9) + G* k (z9)V 2 (z, 9)dz = 0. 

The pair of the equations are defined parameters A k _i(9) and <p k -3(9), where 
k > 2. 

Below we obviously show the process for definition of the parameters for 
A 1 (9) and <p (6). 



17 



The second-order correction term of asymptotic expansion (J31|) . 

The substitution of $^ in the WKB asymptotic expansion leads to the 
equation for the second correction term: 

itfdsQ? +(-9 + 2\U u) \ 2 ^j <# + (u^ 2 - 1^ (Q^)* = G 2 , 

G 2 = - 2\$P\ 2 U® + (<f>?) 2 (UW)* 

The particular solution of the inhomogeneous equation has a form: 

$® = V^S, 9) J (g 2 (ct, 9)V*(a, 9) + G* 2 (a, 9)V 2 (a, 9) 

V 2 (S, 9) J" (g 2 (<t, 9)V*(a, 9) + G*(a, 9)V x (a, 

It is easy to see the solution of the equation is bounded with respect to the 
fast variable S. 



The third-order correction term of asymptotic expansion (13 lj) . The 

equation for the third-order term of the WKB-solution is following 



9 + 2\U {j) \ 2 Wi 



■■Jo) 



G 3 = -id e <f>? ~ ids$?<Po ~ \$?\ 2 ®? ~ 

Equations (|33|) give the differential equations for the parameters Aq(9) 
and fo(9). One can obtain the equations for the main terms of the asymptotic 
of Aq(9) and <fo(9) using the obvious form for the asymptotic expansion of 
U® as e -> 

3A + 4(1 + 9)A' + 0(e) = 0, 
-2(1 + 9) 3/2 ip f + A 2 (16 + 120 - 49 2 + (8 + 89 + 39 2 )V9 + 1) + 0(e) = 0. 

Using regular asymptotic expansion with respect to e we obtain: 



A (9) 



[1 + 9) 3 / 4 ' 



<Po(6) ~ ip 
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Here A 0t0 and <^o,o are parameters of the asymptotic solution for equation (0). 
The values of the parameters will be obtained by matching of the constructed 
WKB-expansion and the asymptotic expansion in the Painleve layer. 

As a result we write out the asymptotic formula for the solution as e — > 0: 

&i\S,9,e) ~ {-iy iv / r + ^ + 5 1/ %,o(^L=cos(S)+z^rT^sin( 1 S)^ 

S ~ j + <p ,o + ^ (59 + + |(112- 86)VT+6 + 3 ln(l + O^J . 
Here A 00 and c^ ,o are the parameters of the constructed WKB-solution. 



The domain of validity for the WKB-asymptotic solution at short 
distance of the slowly varying equilibrium. The constructed WKB- 
asymptotic solution is not valid at the turning points. The higher-order terms 
of the asymptotic expansion are infinite at these points. The turning points 
are 9 = — 1 and 9 = 00. Using obvious formulas for the higher-order terms 
we reduce: 

$0') = 0({9 + l) 3/2 $ Hi)> 0-^-1 + 0, 
Then the domain of validity for the WKB-solution is 

e -2/3 (l + 0) » 1, 9^1 + 0. 

When 9 — > 00 the domain of validity for the WKB-solution is obtained 
analogously: 

*%> = 0(6V'*® 1 ), 9^oo, 

therefore: 

9 < £- 4/5 , 9 -> 00. 



Matching the expansion in the Painleve layer and the captured 
WKB-asymptotic expansion Asymptotic expansion (fT2"|) is valid as 
—e 2 ^ 3 r] < 1 or + 1) C 1 for 9 close to —1. The captured asymptotic 
expansion is valid as e~ 2 ^ 3 (l + 9) ^> 1. Therefore captured asymptotic so- 
lution (|31[) and asymptotic expansion (jl2j) are valid both in the domain 
£ 2 /3 ^ 1 Tj s i n g uniqueness theorem for the asymptotic ex- 

pansions one can obtain that these expansions coincide in this domain. The 
usual matching procedure j20] gives us the connection formulas for the pa- 
rameters A/, and 0^ of expansion ()31|) and parameters of expansion f)12|) . For 
example, the parameters of the first correction term of (|31|) are: 

V? ,o = v, A 0)0 = p. (34) 
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Formulas ([24)1 . (jSlj) and theorem |2] give us the statement of theorem ^ 
Acknowledgments. We are grateful to L. Friedland, R.N. Garifullin, 
L.A. Kalyakin and B.I. Suleimanov for helpful comments and discussions. 
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